Dynamical and photometric imprints of feedback processes on the 
formation and evolution of E/SO galaxies. 
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& ■ ABSTRACT 

We show that the observed Velocity Dispersion Function of E/SO galaxies matches strikingly 
well the distribution function of virial velocities of massive halos virializing at z > 1.5, as predicted 
by the standard hierarchical clustering scenario in a ACDM cosmology, for a constant ratio 
er/Vvir ~ 0.55 ± 0.05, close to the value expected at virialization if it typically occurred at z k> 3. 
This strongly suggests that dissipative processes and later merging events had little impact on 
the matter density profile. Adopting the above cr/Kir ratio, the observed relationships between 
photometric and dynamical properties which define the fundamental plane of elliptical galaxies, 
such as the luminosity-cr (Faber-Jackson) and the luminosity-effective radius relations, as well 
as the Mbh-c relation, are nicely reproduced. Their shapes turn out to be determined by the 
mutual feedback of star- formation (and supernova explosions) and nuclear activity, along the lines 
discussed by Granato et al. (2004). To our knowledge, this is the first semi-analytic model for 
which simultaneous fits of the fundamental plane relations and of the epoch-dependent luminosity 
function of spheroidal galaxies have been presented. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: formation — 
quasars: general 
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Universe from the formation of the first stars in 
the so called "dark ages" to the present epoch, ex- 
plaining the variety of structures from galaxies to 
clusters. Understanding the formation and evo- 
lution of spheroidal galaxies, which comprise the 
most massive and most luminous galaxies in the 
universe and contain a large fraction of all stars, 
turned out to be particularly intriguing. Hierar- 
chical models tend to predict substantially fewer 
massive galaxies at high redshift than are observed 
(Blain et al. 2002; Scott et al. 2002; Daddi 
et al. 2003; Tecza et al. 2004; Somerville et 
al. 2004). On the other hand, the traditional 
monolithic models, whereby these objects formed 
their stars on a timescale shorter than their free- 
fall time and evolved passively thereafter (Eggcn, 
Lynden-Bell & Sandage 1962), do not fit into a co- 
herent scenario for structure formation from pri- 
mordial density perturbations, and also tend to 
over-predict high-redshift galaxies. 

In this situation it is essential to look for guid- 
ance from observational data. It has long been 
known that stellar populations in elliptical galax- 
ies are old and essentially coeval (Sandage & 
Visvanathan 1978; Bernardi et al 1998; Trager 
et al. 2000; Terlevich & Forbes 2002). A 
color-magnitude relation is also well established: 
brighter spheroids are redder (Bower et al. 1992). 
The widely accepted interpretation is that brighter 
objects are richer in metals and the spread of their 
star formation epochs is small enough to avoid 
smearing of their colors. The slope of this relation 
does not change with redshift (Ellis et al. 1997; 
Kodama et al. 1998) supporting this interpre- 
tation. The star formation history of spheroidal 
galaxies is mirrored in the Fundamental Plane 
(Djorgovski & Davies 1987; Dressier et al. 1987) 
and in its evolution with redshift. Elliptical galax- 
ies adhere to this plane with a surprisingly low 
orthogonal scatter (~ 15%), as expected for a ho- 
mogeneous family of galaxies. Recent studies (e.g. 
Treu et al. 2002; van der Wei et al. 2004; Holden 
et al. 2004, 2005) suggest that ellipticals, both in 
the field and in clusters, follow this fundamental 
relation up to z ~ 1, consistent with the hypothe- 
sis that massive spheroids are old and quiescent. 

Direct evidence that massive galaxies with M > 
10 n M Q were in place at z > 2, is provided 
by recent if -band surveys (Cimatti et al. 2002; 
Kashikawa et al. 2003; Fontana et al. 2005). The 



space density of Extremely Red Objects (EROs) 
at z > 3 is only a factor ~ 5-10 less than at z <~ 1 
(Tecza et al. 2004). 

While the evolution of dark matter halos is con- 
trolled only by gravity, and therefore the underly- 
ing physics is simple (even if the evolutionary be- 
haviour is complex), processes involving baryons 
are intricate and may hold the key to reconcile 
theory with observations. As shown by Granato 
et al. (2004) the feedback from supernova (SN) ex- 
plosions and from nuclear activity can reverse the 
hierarchical scenario for baryons (see also Granato 
2001). In other words, the canonical hierarchi- 
cal Cold Dark Matter (CDM) scheme - small 
clumps collapse first - is reversed for baryon col- 
lapse and the formation of luminous objects (Anti- 
hierarchical Baryon Collapse (ABC) scenario). 

In this paper, after a short overview of the 
Granato et al. (2004) model (Sect. 2), we in- 
vestigate the imprints of processes governing the 
formation and the early evolution of spheroidal 
galaxies on their local dynamic and photometric 
properties, such as the Velocity Dispersion Func- 
tion (VDF, Sect. 3.1), the luminosity- velocity dis- 
persion relation (Faber & Jackson 1976, Sect. 3.2), 
the stellar luminosity-effective radius relation 
(Bernardi et al. 2003, Sect. 3.3), and the BH mass- 
velocity dispersion relation (Ferrarese & Merritt 
2000; Tremaine et al. 2002; Onken et al. 2004; 
Sect. 3.4). The main conclusions are presented 
and discussed in Sect. 4. 

Throughout this work we adopt a spatially 
flat cold dark matter cosmology with cosmolog- 
ical constant, consistent with the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) data (Ben- 
nett et al. 2003): tt m = 0.3, tt b = 0.047, and 
tt A = 0.70, H = 70 km s" 1 Mpc"\ cr 8 = 0.84, 
and an index n = 1.0 for the power spectrum of 
primordial density fluctuations. 

2. Overview of the Granato et al. (2004) 
model 

While referring to the Granato et al. (2004) 
paper for a full account of the model assumptions 
and their physical justification, we provide here, 
for the reader's convenience, a brief summary of 
its main features. 

The model follows with simple, physically 
grounded recipes and a semi-analytic technique 
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the evolution of the baryonic component of proto- 
sphcroidal galaxies within massive dark matter 
(DM) halos forming at the rate predicted by the 
standard hierarchical clustering scenario within a 
ACDM cosmology. The main novelty with respect 
to other semi-analytic models is the central role 
attributed to the mutual feedback between star 
formation and growth of a super massive black 
hole (SMBH) in the galaxy center. 

The idea that SN and QSO feedback play an im- 
portant role in the evolution of spheroidal galaxies 
has been pointed out by several authors (Dekel 
& Silk 1986; White & Frenk 1991; Haehnelt, 
Natarajan & Rees 1998; Silk & Rees 1998; Fabian 
1999). Granato et al. (2004) worked out, for 
the first time, the symbiotic evolution of the host 
galaxy and the central black-hole (BH), includ- 
ing the feedback. In this model, the formation 
rate of massive halos (2.5 x 10 n M Q < M vir < 
2 x 10 13 M Q ) is approximated by the positive part 
of the time derivative of the halo mass function 
(Press & Schechter 1974, revised by Sheth & Tor- 
men 2002). The gas, heated at virial temperature 
and moderately clumpy (clumping factor ~ 20), 
cools to form stars, especially in the innermost re- 
gions where the density is the highest. The radia- 
tion drag due to starlight acts on the cold gas, fur- 
ther decreasing its angular momentum and caus- 
ing an inflow into a reservoir around the central 
BH, to be subsequently accreted into it, increas- 
ing its mass and powering the nuclear activity. 
In turn, the feedbacks from SN explosions and 
from the active nucleus regulate the star forma- 
tion rate and the gas inflow, and eventually un- 
bind the residual gas, thus halting both the star 
formation and the BH growth. In fact, important 
parameters are the efficiency of SN energy transfer 
to the cold gas (csn) and the fraction of the QSO 
luminosity in winds. 

Further relevant differences with other semi- 
analytic models are the allowance for a clumping 
factor, C ~ 20, speeding up the radiative cool- 
ing so that even in very massive halos (M v - U ~ 
10 13 M Q ), the gas, heated to the virial tempera- 
ture, can cool on a relatively short (~ 0.5-1 Gyr) 
timescale, at least in the dense central regions, 
and the assumption that the large-scale angular 
momentum does not effectively slow down the col- 
lapse of the gas and the star formation in massive 
halos virializcd at high rcdshift. As for the latter 



assumption, we note that during the fast collapse 
phase, when the potential well associated with a 
galactic halo is established, gas clouds can lose 
their orbital energy to the dark matter by dynam- 
ical friction on a timescale shorter than the halo 
collapse timescale (see, e.g., Mo & Mao 2004). 

The model prescriptions are assumed to ap- 
ply to DM halos virializing at z v ; r 1.5 and 

M vir > 2.5 x lO n M0. These cuts are meant to 

crudely single out galactic halos associated with 
spheroidal galaxies. Disk (and irregular) galax- 
ies are envisaged as associated primarily to halos 
virializing at z v - K £ 1.5, some of which have incor- 
porated most halos less massive than 2.5 x 10 11 M 
virializing at earlier times, that may become the 
bulges of late type galaxies. However, the model 
does not address the formation of disk (and ir- 
regular) galaxies and these objects will not be 
considered in this paper. 

The kinetic energy fed by supernovae is in- 
creasingly effective, with decreasing halo mass, in 
slowing down (and eventually halting) both the 
star formation and the gas accretion onto the cen- 
tral black hole. On the contrary, star formation 
and black hole growth proceed very effectively in 
the more massive halos, giving rise to the bright 
SCUBA phase, until the energy injected by the 
active nucleus in the surrounding interstellar gas 
unbinds it, thus halting both the star formation 
and the black hole growth (and establishing the 
observed relationship between black hole mass and 
stellar velocity dispersion or halo mass). Not only 
the black hole growth is faster in more massive ha- 
los, but also the feedback of the active nucleus on 
the interstellar medium is stronger, to the effect 
of sweeping out such medium earlier, thus caus- 
ing a shorter duration of the active star-formation 
phase. 

The basic yields of the model are the star- 
formation rate, ip(t), as a function of the galactic 
age, t, (hence the evolution of the mass in stars, 
M S ph(t)), and the growth of the central BH mass, 
MBn(t), for any given value of the halo mass, M v ; r , 
and of the virialization rcdshift, z v ; r . These quan- 
tities are obtained solving the system of differen- 
tial equations given by Granato et al. (2004), i.e. 
their eqs. (10), (16), and (23). The rates at which 
the diffuse gas cools [M C oid(i)] and forms stars are 
given by their eqs. (9) and (10), respectively, while 
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their eq. (16) yields the growth rate - due to 
the radiation drag mechanism - of the central BH 
reservoir [-Minflow(i)], allowing the building-up of 
the BH mass (eq. 23). Finally, the rates at which 
the feedback from SN [M c s Y d (t)] and QSO [M?f°] 
heats the cold gas moving it into the hot phase are 
given by their eqs. (11) and (31), respectively. In 
their notation: 

Minf(t) = -Mcoi d (t) + M%f° (1) 

M gas (t) = (t) + M cold , sf (t) (2) 

M ICS (t) = M inflow (t) - M B h(*) (3) 
Afcoid, rf (t) = Mcoid(t) - V>(0 + M™ ld (t) 

+M?™(t), (4) 

where a dot above a symbol denotes a time 
derivative. Note that, in eq. (10) of Granato et 
al. (2004), defining ip(t), M co id(r, i) is actually 

Mcold.sf (»",*). 

Once the star- formation history, ip(t), of a 
galaxy has been computed, its luminosity in 
any chosen band is obtained, as a function of 
t, from the GRASIL code (Silva et al. 1998), 
which, given yields the chemical and the 

spectro-photometric evolution from the radio to 
the X-ray band, allowing for the effect of dust 
absorption and reradiation. GRASIL is avail- 
able at web.pd.astro.it/granato and at adlibi- 
tum.oat.ts.astro.it/silva/default.html. As men- 
tioned above, the halo formation rate is obtained 
using the Press & Schechter (1974) formalism, as 
improved by Shcth & Tormen (2002; see eq. (5) of 
Granato et al. 2004). 

Putting these ingredients together we obtain 
the model predictions for the statistical proper- 
ties of spheroidal galaxies as a function of cosmic 
time. Although the model has several parame- 
ters (listed in Table 1 of Granato et al. 2004), 
all of them are strongly constrained by the wealth 
of observational data that have become available 
in recent years, including the luminosity functions 
of spheroidal galaxies at different redshifts, their 
counts and rcdshift distributions in optical and 
near-IR bands, the 850 ^m SCUBA counts and the 
associated (if preliminary) redshift distributions, 
the relationship between the black-hole mass and 
the stellar velocity dispersion, the local black-hole 
mass function, and more (Granato et al. 2004; 
Silva et al. 2005; Shankar et al. 2004). 
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Fig. 1. — Comparison of the virial velocity func- 
tion of galaxies with z v - u > 1.5 implied by the stan- 
dard gravitational clustering scenario (solid line, 
upper scale), with observational estimates of the 
velocity dispersion function of early-type galaxies 
by Shankar et al. (2004; shaded area) and Sheth 
et al. (2003; triple-dot dashed line). The crosses 
show the contribution of bulges of late- type galax- 
ies, determined by Sheth et al. (2003). The two 
functions match for a = 0.55Vv; r . The dashed 
lines show the virial velocity functions of galaxies 
for different choices of the minimum z v i r : from top 
tO bottom, Z v ir,min = 2, 3, 4, 5. 

3. Properties of early type galaxies 
3.1. The Velocity Dispersion Function 

As discussed by Loeb & Peebles (2003), the lo- 
cal Velocity Dispersion Function (VDF) can pro- 
vide interesting hints on the structure formation 
process 1 . Accurate determinations of the VDF of 
early-type galaxies have been obtained by Shcth 
et al. (2003) and Shankar et al. (2004), based 
on a large sample (~ 9000 E/S0 galaxies) drawn 
from the SDSS (Bcrnardi et al. 2003). The two 
determinations agree remarkably well (see Fig. 1) 
except for low values of the velocity dispersion, 
where the Sheth et al. (2003) VDF declines due 
to the selection criteria adopted by Bernardi et al. 
(2003) to compute the low luminosity portion of 
the luminosity function. 

As for dark matter halos, it is convenient to 

Throughout this paper the line-of-sight velocity dispersion 
a is referred to R e /8, where the effective radius R e is de- 
fined as the radius containing half of the light. 
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define a 'viriaP velocity, equal to the circular ve- 
locity at the virial radius (Navarro et al. 1997, 
NFW; Bullock et al. 2001): 

Kl = ^ • (5) 
Since, given the virialization redshift, Kir depends 

1 /3 

only on M v i r (Kir oc Af ir ) the Kir distribution 
function can be straightforwardly derived from 
the mass distribution function of spheroidal galax- 
ies [eq. (5) of Granato et al. (2004)], integrated 
over the virialization redshifts. Following Granato 
et al. (2004), we assume that all massive halos 
(2.5 x 10 n M Q < M vil < 2 x 1O 13 M ) virializ- 
ing at z > 1.5 yield spheroidal galaxies or bulges 
of later type galaxies. It may be noted, in pass- 
ing, that the adopted upper mass limit is close to 
that inferred by Kochanek & White (2001) from 
the distribution of gravitational lens image sepa- 
rations. 

As illustrated by Fig. 1, the derived Kir distri- 
bution function (which depends only on the evo- 
lution of dark matter halos) accurately matches 
the observationally determined VDF (which may 
be affected by the physics of baryons, and in par- 
ticular by dissipative processes) if a ~ 0.55Kir- 
The best fit and the confidence intervals of the 
c/Kir ra tio depend somewhat on the choice of the 
upper mass limit, M sup (reference value M sup = 
2 x 10 13 M o ) and of cr 8 (reference value cr 8 = 0.84), 
while the choice of the minimum virialization red- 
shift (z v ir,min) set at 1.5, does not affect apprecia- 
bly the fit. For cr 8 = 0.84, the 95% confidence 
interval, determined utilizing the x 2 statistic with 
2 degrees of freedom, is 0.53 < a /Kir < 0.60, 
the upper end corresponding to the smallest value 
of M sup , which is constrained to be > 1O 13 M 
to ensure consistency with the local if -band lu- 
minosity function (see Granato et al. 2004). De- 
creasing as to 0.8, yields a 95% confidence interval 
0.50 < a/Vvir < 0.61, while increasing erg to 0.9, 
we get an acceptable fit for M sup ~ 10 13 M Q or 
cr/Kir — 0.50. As already mentioned, the data 
indicate a a /Kir ratio independent of virial mass. 
A weak dependence, not steeper than cx/Kir °c 
M°; r 05 , is however allowed. 

A linear relationship between the central veloc- 
ity dispersion, cr, and the maximum circular veloc- 
ity, w ™ ax - which, for a concentration c ~ 3 (see 
below), is essentially equal to Kir _ was reported 



by Gerhard et al. (2001) for a sample of 21 mostly 
luminous, slowly rotating elliptical galaxies, al- 
though the ratio is somewhat higher than found 
here: a = 0.66t>™ ax . A weakly non-linear rela- 
tionship was found by Ferrarese (2002) for a sam- 
ple of 13 spiral galaxies with rotation curves ex- 
tending beyond the B = 25 mag/arcsec 2 isophote: 
cr/200kms- 1 = O.eO^^OOkms" 1 ) 119 . Ac- 
cording to Ferrarese (2002) this relation can be 
considered valid also in the a range populated by 
elliptical galaxies. 

If the stellar velocity dispersion profile is ap- 
proximately isothermal and stellar velocities are 
isotropic, adopting the Navarro et al. (1997) den- 
sity profile we obtain the following relationship be- 
tween Kir and the velocity dispersion a: 

a _ {c[3c 2 + 4c - 2cln(l + c) - 21n(l + c)]} 1 / 2 
V~'"' V6[l + (l + c) ln(l + c)] 

(6) 

where c, equal to the ratio of i? v ; r to the NFW in- 
ner radius r s , is the "concentration" . The N-body 
simulations by Zhao et al. (2003a) show that halos 
of mass greater than 10 11 h _1 M Q at z > 3 have all 
a similar median concentration c ~ 3.5. For c = 2, 
3, 4, eq. (6) yields a/Kir = 0.49, 0.57, 0.62, re- 
spectively. Thus, the value of cr/Kir for which we 
get a match between the local VDF and the Kir 
distribution function of dark halos is remarkably 
close to the value expected based on simulations. 

The tight correspondence between the VDF 
and the velocity distribution function of dark ha- 
los lends support to the dynamical attractor hy- 
pothesis (Loeb & Peebles 2003; Gao et al. 2004), 
according to which the total distribution of colli- 
sionless matter (dark matter plus stars) keeps es- 
sentially constant in the presence of merging and 
of dissipative settling of baryons, with the dark 
matter distribution expanding to compensate for 
the dissipative settling of baryons (but see Gnedin 
et al. 2004 for a different conclusion). 

The stability of the stellar dynamics in the cen- 
tral regions of dark halos against merging events 
subsequent to the virialization redshift is also con- 
sistent with the results of detailed numerical sim- 
ulations (Wechsler et al. 2002; Zhao et al. 2003b) 
showing that the halo circular velocity changes 
very little after the end of the initial fast accretion 
process, during which most of the specific binding 
energy is assembled, even though a large fraction 
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Fig. 2. — Observed Faber- Jackson relation (filled 
dots) from Bernardi et al. (2003), compared with 
our predictions for different virialization rcdshifts 
z v i r . The curves at the bottom of the figure rep- 
resent the normalized distributions of velocity dis- 
persions of galaxies in 4 absolute magnitude bins 
0.5 mag wide, centered at Mr* = —20.2, —21, 
—22, and —23 (from left to right), as predicted by 
the ABC model. The FWHMs of the predicted 
distributions are remarkably close to the observed 
values (FWHM - 0.09; Bernardi et al. 2003). 

of the halo mass is acquired during the subsequent 
prolonged slow accretion phase. 

Figure 1 also shows the contributions to the 
VDF of different virialization rcdshifts and high- 
lights that the highest velocity portion comes from 
the highest virialization redshifts. This result is 
nicely consistent with the findings by Loeb & Pee- 
bles (2003) who computed the expected cumula- 
tive comoving VDF at z = 4 for the NFW and 
the Moore et al. (1999) density profile, and found 
it consistent with the observational determination 
by Shcth et al. (2003) for velocity dispersions 
a > 300 km s~ 1 , while for lower values of a the pre- 
dicted VDF is low compared to the observed one. 
Our curve for z v ; r > 4 has a very similar behavior; 
but when we allow also for spheroidal galaxies or 
bulges virializing at lower redshifts we can fully ac- 
count for the observed VDF. As mentioned above, 
the contributions of sources virializing at z < 1.5 
is negligibly small. 



3.2. The Faber - Jackson relation 

Since the first measurements of velocity disper- 
sions of early type galaxies were made, it was rec- 
ognized that they are correlated with the galaxy 
luminosities (Poveda 1961; Minkowski 1962). 
Faber & Jackson (1976) showed that Lb oc cr 4 
(Faber- Jackson relation). Bernardi et al. (2003), 
using a sample of ~ 9000 early type galaxies drawn 
from the Sloan Digital Sky Survey (SDSS) in the 
redshift range 0.01 < z < 0.3, found L r , cx cr 3 92 , 
consistent with previous studies (Forbes & Pon- 
man 1999; Pahre et al. 1998). Their data in u, g, 
i, and z bands show that the relation is roughly 
independent of wavelength. The distributions of 
a at fixed luminosity arc approximately Gaussian. 

A detailed quantitative analysis can be carried 
out using the ABC model to follow the time evo- 
lution of the baryonic component, both in the 
gas and in the stellar phase, of each halo mass 
M v i r since the virialization redshift z v - lr . Then, 
the present-day luminosity in any chosen band 
bands from X-ray to radio can be computed with 
the spectrophotometric code GRASIL (Silva et al. 
1998) as a function of V v i r , which is uniquely de- 
termined by the value of M v j r , given z v j r . As il- 
lustrated by Fig. 2, the model not only predicts 
the correct slope of the Faber-Jackson relation, 
but also the correct normalization, when we use 
the ratio a/V v -„ = 0.55 found from the analysis 
of the VDF. Acceptable fits can be obtained with 
0.50 < (T/V v i r < 0.65 (again, the confidence inter- 
val is derived using the \ 2 statistic). 

The scatter in the observed relation is inter- 
preted as an intrinsic property of elliptical galax- 
ies, accounted for by different virialization epochs: 
galaxies with the same spheroidal luminosity, but 
virializing at lower redshifts, have lower velocity 
dispersions. The curves at the bottom of Fig. 2 
show the cr distributions (arbitrary units) of galax- 
ies in four luminosity bins, as obtained from the 
Granato et al. (2004) model. Such distributions 
have been computed integrating over cosmic time, 
for each value of a, the formation rate of halos 
with present day luminosities within the consid- 
ered bin; they turn out to be roughly Gaussian, 
with the peak close to the a expected from the 
best fit relation and FWHM in agreement with 
that observed by Bernardi et al. (2003). 

It is worth noting that the standard scaling of 
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the virial parameters in the hierarchical cluster- 
ing scenario gives (Bullock et al. 2001) M v ; r oc 
V^ r (l + z v i r )~ 3 / 2 , which, for a roughly constant 
M v j r /L ratio, would imply a flatter slope than 
is observed in the Faber-Jackson relation. How- 
ever, the slope is steepened in the ABC model 
which predicts a decrease of the M v i r /M spri ra- 
tio (M sp h being the mass in stars), with increas- 
ing Mvi r , whose details depend on the virializa- 
tion redshift (see Fig. 5 of Granato et al. 2004). 
This is due to feedback from supernovae, which is 
increasingly efficient with decreasing M v ; r in pre- 
venting the gas from cooling and forming stars, 
tempered by the feedback from active nuclei which 
is more effective in the more massive objects. On 
average, we have, to a sufficient approximation, 
M v i r /M S p h oc M^ /5 . The ABC model also pre- 
dicts an essentially constant M spri /£ ratio (the ob- 
served weak luminosity dependence is attributed 
to the systematic changes with luminosity of the 
galactic structure, see Sect. 3.3). Therefore, the 
M v i r -Vvi r relation, for cr/Vvir = const, translates 
into L oc cr 18 / 5 (l + z vir )~ 9 / 4 . The mean slope 
steepens somewhat, approaching that of the ob- 
served Faber-Jackson relation, when we account 
for the varying contributions of different virial- 
ization rcdshifts to different a intervals. A fur- 
ther steepening of the Faber-Jackson relation is 
expected at low a values, corresponding to less 
massive objects where the SN feedback, yielding 
Afvir/Mjph oc M vi ^ 2 , dominates. 

In conclusion, the Faber-Jackson relation is in- 
terpreted as providing a quantitative measure of 
the effect of feedback, and primarily of the effect 
of the energy injected onto the interstellar medium 
by supernovae. The close agreement with the pre- 
dictions of the ABC confirms the correctness of 
the adopted recipes. We note that fitting simul- 
taneously the galaxy epoch dependent luminosity 
functions and dynamical properties has long been 
a very challenging problem for semi-analytic mod- 
els (see, e.g., Mo & Mao 2004). In the case of 
spheroidal galaxies, this requires that the model 
correctly predicts not only the evolution oiL/M vlr 
but also the distribution of virialization rcdshifts 
as a function of the present-day luminosity. To our 
knowledge, this is the first semi-analytic model for 
which successful fits of all these quantities have 
been reported. 
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Fig. 3. — Absolute magnitude in the r* band ver- 
sus R e . The data points are from Bernardi et al. 
(2003) while the lines show our model predictions 
for different virialization rcdshifts. 

3.3. Completing the Fundamental Plane 

An additional, though not independent, check 
on the physical processes involving baryons in the 
ABC model is provided by the observed L-R e rela- 
tion (see e.g. Bernardi et al. 2003). The effective 
radius, like the velocity dispersion, is related to 
the collapse and settling of the baryonic compo- 
nent inside the DM potential well. On the other 
hand, the model does not give us dynamical infor- 
mation. However, it yields the total mass in stars 
and dark remnants, while observations ensure that 
the starlight distribution in spheroids has quite a 
uniform shape. In fact the surface brightness dis- 
tribution is well represented by log(/(r)) oc r 1 /" 
(Sersic 1968); the classical de Vaucouleurs profile 
has n — 4 (de Vaucouleurs 1948). 

Bordello et al. (2003), properly taking into ac- 
count the light distribution, the mass traced by 
light, and the dark matter, found that the effec- 
tive radius can be written as: 

fl e = (fc ff +a DM )^* (7) 
a 

where M spri is the total mass traced by light, k a 
is a constant depending on the light distribution 
(fc CT =0.174 for the de Vaucouleurs profile) and cudm 
is a function of the DM mass inside R e and of 
its density distribution [see eqs. (13) and (14) of 
Borricllo et al. 2003]. As found by Gerhard et al. 
(2001) for a sample of about 20 elliptical galaxies 
and generalized by Bordello et al. (2003) on the 
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basis of the narrowness of the Fundamental Plane, 
the DM inside the effective radius amounts to only 
10-50% of the total mass. Such an amount and 
its distribution make «dm almost negligible in the 
above equation, that can be rearranged as follows, 
using the definition of V V - 1T [eq. (5)]: 
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Re 
Rv'u 



(k a + qdm) 



M sph /M v 
(a/Kir) 1 



(8) 



The ABC model yields the M sp h/M v i r ratio and 
the ensuing M sp h/L ratio, whose value depends on 
the adopted IMF. Using the Salpeter's IMF we get 
M sp h/L r ~ 5M q /Lq, with very little dependence 
on the luminosity in the reference r* band. 

On the other hand, systematic changes with lu- 
minosity of the galactic structure, quantified by 
the Scrsic index n, have been reported and found 
to account for the variation of the M/L ratio with 
luminosity (Ciotti et al 1996; Graham et al. 2001; 
Trujillo et al. 2004). Bcrnardi et al. (2003), 
exploiting the large sample of early-type galax- 
ies in the SDSS, defined an effective mass M = 
4i? e cr 2 /G, and found M /L r = 3.5(L r /L*) 015 , in 
solar units, with L* = 2 x 10 10 L Q . After eq. (7) 
M = 4(fc CT + a D M)M sp h, so that, neglecting a DM 
and setting M sp h/L r — 5M Q /L Q , we have: 



fc CT ~0.174(-^) 



r NO. 15 



(9) 



consistent with the findings of Borriello et al. 
(2003) based on a much smaller sample (221 
nearby galaxies) but with more detailed obser- 
vations. Hence: 



^ „ „_ ,L r . 15 L r /M v j T 

Rvir " ' y L*' (a/Kir) 2 ' 



(10) 



The ABC model gives, for each value of z v ir and 
for an assumed IMF, i? v ir and L r function 
of M v ir- Using the above equation we can then 
obtain a relationship between L r (or the abso- 
lute magnitude M r *) and R e . The results for a 
Salpeter IMF and a/Kir = 0.55, as implied by 
the VDF, are compared in Fig. 3 with the data of 
Bernardi et al. (2003) . The presence of an average 
-30% of DM inside R e would imply a DM < 0.02 
(see Borriello et al. 2003) and would not modify 
the fit. 

The M r *-R e relation is related to the Faber- 
Jackson relation through eq. (7), given the ob- 
servationally determined L r /M sp h ratio. Thus, 
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Fig. 4. — The Mbh — cr relation predicted by 
the model for different virialization redshifts com- 
pared with observational data. Filled circles rep- 
resent nearby galaxies with dynamical measure- 
ments of the BH mass (Tremaine et al. 2002) and 
the upper limits for M33 (Gebhardt et al. 2001) 
and NGC 205 (Valluri et al. 2005). Open symbols 
refer to galaxies with BH masses estimated from 
the H/3 linewidth-luminosity-mass relation: NGC 
4395 (open diamond; Filippenko & Ho 2003), POX 
52 (open circle; Barth et al. 2004) and a sample of 
7 faint active nuclei drawn from the SDSS, (open 
squares; Greene et al. 2004). 



it does not provide an independent test of the 
L r /M v i r ratio predicted by the ABC model. Still, 
the consistency of the model with a different set 
of data is, at least, reassuring. 

3.4. The central black hole 

In the last years massive BHs were proven to 
be ubiquitous in massive E/S0 galaxies (see e.g. 
Magorrian et al. 1998). Tight relations have been 
discovered between BH masses, Mbh, and stel- 
lar velocity dispersions (Ferrarese & Merritt 2000; 
Gebhardt et al. 2000; Tremaine et al. 2002; Onkcn 
et al. 2004) , masses of the spheroidal components 
(McLure & Dunlop 2002; Dunlop et al. 2003; Mar- 
coni & Hunt 2003), and masses of the dark halos 
(Ferrarese 2002); see Ferrarese & Ford (2004) for 
a comprehensive review. The Mbh _ ct relationship 
has been extensively discussed as a probe of the 
interplay between QSO and host galaxy evolution 
(Silk & Rees 1998; Fabian 1999; Cavaliere et al. 
2002). 
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Fig. 5. Predicted Mbh -M sp h relation for dif- 
ferent virialization rcdshifts compared with obser- 
vational data from Hiring & Rix (2004). Open 
squares are sources from their Group 1 and solid 
circles from their Group 2. 

Granato et al. (2004) showed that the growth 
of the stellar component and of the central BH 
are strictly symbiotic, double tied by the posi- 
tive effect of the photon drag, proportional to the 
star formation rate, favoring the inflow toward the 
reservoir around the central BH (Umemura 2001; 
Kawakatu & Umemura 2002), and by the negative 
effect of the BH growth, powering outflows. 

The ABC model follows, for any given M v i r and 
z v i r , the growth of the central BH and gives its 
final mass. The predicted M B H _ Kir relation, at 
fixed z v ; r , is immediately translated into Mbh-c 
using the ratio cr/Kir = 0.55 derived from the 
VDF. The scatter of data points is interpreted as 
reflecting the distribution of z v - lr . 

The model also predicts a steepening of the re- 
lation at low a values, due to the combined ef- 
fect of SN feedback - which is increasingly effi- 
cient with decreasing halo mass in slowing down 
the gas infall onto the central BH - and of the de- 
creased radiation drag (due to a decrease of the 
optical depth). From an observational point of 
view, the behavior of the Mbh— c relation in the 
low BH mass and low velocity dispersion regime 
is still unclear, due to the dearth of data and 
to the uncertainties on the Mbh estimates. In 
Fig. 4 we plot objects found in literature with 
estimated BH masses Mbh ^ 10 6 M Q . Dynami- 
cal measurements are available only for M33 and 



NGC 205, while all the other BH masses have 
been estimated through the linewidth-luminosity- 
mass scaling relation (Kaspi et al. 2000). There- 
fore the cases of M33 and NGC 205 are partic- 
ularly interesting, although they are outside the 
range of masses to which we applied our model 
up to now. For M33 the upper limit to the BH 
mass is M BH ^ 3000M Q (Merritt et al. 2001) or 
M BH ^ 1500M Q (Gebhardt et al. 2001), while, for 
NGC 205, M BH ^ 2.2 x 10 4 M Q . Both fall below 
the extrapolation of the M^n-a relation holding 
at higher a (Ferrarese & Merritt 2000; Gebhardt 
et al. 2000; Tremaine et al. 2002; Onkcn et al. 
2004), while it is fully consistent with the extrap- 
olation of our model. On the other hand, the BH 
mass estimates for the faint type 1 Seyfert nuclei 
in NGC 4395 (Filippenko & Ho 2003), in POX 52 
(Barth et al. 2004) and in 7 galaxies drawn from 
SDSS (Greene et al. 2004) are only marginally 
consistent with the steepening predicted by our 
model. More data on the so called intermedi- 
ate mass BHs (1O 3 M < M BH < 10 6 M Q ) in the 
galactic centers (see van der Marel (2003) for a re- 
view) are needed to clarify this issue. Of course 
we should also keep in mind that the low BH mass 
portion of the diagram might just reflect the distri- 
bution of BH seeds (possibly created by merging 
of smaller BHs during the fast accretion epoch) 
and be only weakly affected by the mass accre- 
tion which is controlled by the effects discussed 
above. For example, Koushiappas et al. (2004) 
presented a model yielding seed BHs with char- 
acteristic masses ~ 10 5 M . Also, low BH masses 
may be easily increased by substantial factors dur- 
ing later accretion phases, when host galaxies are 
disk- (rather than bulge-) dominated. 

Another prediction of the ABC model is that 
the ratio of the BH mass to the mass in stars, 
M B H/M sp h, is almost insensitive to variations of 
z v i r . Therefore the M B H _ M sp } 1 relation is expected 
to have a smaller intrinsic scatter (see Fig. 5) than 
the Mbh-ct relation. This is because, as men- 
tioned above, the growth of the BH mass is con- 
trolled by the star formation rate through the ra- 
diation drag and the SN feedback, and in turn, 
the feedback from the active nucleus can eventu- 
ally sweep out the gas thus halting both the star 
formation and the accretion on the BH. Thus the 
stellar and BH mass grow (and stop growing) in 
parallel. The parallelism is not exact, however, 
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since the star formation rate has a twofold effect 
on the BH growth. As a result, the MBH^M sp h 
relation is not strictly linear, but bends down at 
small masses, and is slightly different for differ- 
ent values of z v ; r . On the other hand, estimates 
of M sp h are somewhat indirect and therefore li- 
able to larger uncertainties than those of a, that 
can be directly measured; this may translate in an 
observed scatter around the mean Mbh-MspIi re- 
lation comparable to, or larger than that for the 
Mbh-c relation, in spite of the smaller intrinsic 
scatter. 

In Fig. 5 we plot the Mbh -M sp h relation for 
objects with reliable bulge mass determinations 
from Haring & Rix (2004). The agreement be- 
tween the data, suggesting a slightly non-linear 
relation (Mbh oc M^ 2 ), with our predictions is 
remarkably good. As pointed out by Haring & Rix 
(2004) , a significant fraction of the scatter (which 
is < 0.3 dex) can be attributed to measurement 
errors. Marconi & Hunt (2003) found that the 
scatter in this relation is reduced to ~ 0.25 dex, 
when Mgph is estimated MS c\ virial mass (~ R e o~ 2 ). 

The very good fits of the observed Mbh~c and 
Mbh -M sp h relations are additional strong indica- 
tions that the ABC model properly includes the 
mutual feedbacks of stars and QSOs. It is worth 
noticing that the model also correctly predicts 
the local BH mass function in spheroidal galax- 
ies (Granato et al. 2004; Shankar et al. 2004). 

4. Discussion and conclusions 

We have pointed out an impressively, and unex- 
pectedly, tight correspondence between the virial 
velocities, Kdr, controlled by the dynamics of 
dark halos, and the stellar velocity dispersions, 
a, that feel the effect of dissipative baryon load- 
ing. A straightforward comparison of the virial 
velocity distribution, implied by the standard hi- 
erarchical clustering scenario in a ACDM cosmol- 
ogy, with the observed velocity dispersion function 
of spheroidal galaxies, shows that the two func- 
tions match remarkably well for a constant ra- 
tio cr/Kir — 0.55. For a NFW density profile, 
this ratio corresponds to a concentration param- 
eter c ~ 3, noticeably close to that found from 
N-body simulations (Zhao et al. 2003a) to apply 
to halos of mass greater than 10 11 h _1 M Q viri- 
alizing at z > 3, with no significant dependence 



on halo mass. Thus the matter density profile at 
virialization appears to be essentially unaffected 
by the subsequent events, including mergers and 
the dissipative contraction of baryons, consistent 
with the dynamical attractor hypothesis (Loeb & 
Peebles 2003; Gao et al. 2004). 

The Kir-c relation is a key ingredient to con- 
nect theoretical predictions with observations. Us- 
ing the above determination of the cr/Kir ra- 
tio, we have shown that the observed relation- 
ships between photometric and dynamical prop- 
erties of spheroidal galaxies, defining the Funda- 
mental Plane, carry clear imprints of the feedback 
processes ruling the early evolution of spheroidal 
galaxies. The steeper slope of the luminosity- 
a (Faber-Jackson) relation compared to the pre- 
dicted M v i r -y v i r relation is interpreted as due pri- 
marily to the heating of the interstellar medium 
by supernovae, which increasingly hampers star- 
formation in smaller and smaller halos, and to 
winds driven by the active nucleus, that eventu- 
ally sweep out the residual interstellar gas in the 
most massive nuclei. A full analysis shows that the 
treatment of feedback adopted in the ABC model 
nicely accounts not only for the slope, but also for 
the normalization of the Faber-Jackson relation 
(see Fig. 2). The model implies that the observed 
scatter is mostly intrinsic and due to the spread 
of virialization rcdshifts. It also predicts a mod- 
erate steepening of the relation at low masses. To 
our knowledge, this is the first semi-analytic model 
for which simultaneous fits of the epoch dependent 
luminosity functions and of the Faber-Jackson re- 
lation have been reported. 

The effect of feedback also determines the slope 
of the luminosity-effective radius relation. As dis- 
cussed in Sect. 3.2, the feedback yields, approx- 

5/6 

imately, M V - 1T cx M sph . Then, taking into ac- 
count eq. (8), neglecting ojdm and with k a cx 
L°- 15 , we get, from M vir oc i?^ ir (for given z vir ), 

i? e cx L®: 15 R®(^ and L r cx i? 2 . Weighting with the 
redshift distributions appropriate for each value 
of R e (high-z galaxies contribute more to low 
Re's), the mean L r -R e relationship flattens to- 
wards the observed relation L r cx R^ 58 (Bernardi 
et al. 2003). A further flattening is predicted at 
low R e values, where the SN feedback dominates, 

2/3 

yielding, approximately, M v -„ cx M sphl whence 
R e cx L°- 15 ilV 2 and L r cx R\ A . Again, the model 
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Fig. 6. — Effect of feedback on the Faber- Jackson relation (left panel column), on the local r*-band luminosity 
function of spheroidal galaxies (central panel column), and on the Mbh _ c relation (right panel column), 
according to the ABC model. The upper, central, and lower panel rows compare with the data the ABC 
model predictions with feedback fully included and switching off, in turn, the SN and the QSO feedback, 
respectively. Symbols are as in Figs. 2 and 4. The observational determinations of the local luminosity 
function are from Nakamura et al. (2002; open squares), and Bernardi et al. (2003; filled circles). 



correctly reproduces not only the slope but also 
the normalization of the relation (see Fig. 3). As in 
the case of the Faber- Jackson relation, we expect 
a scatter comparable to the observed one due to 
the different virialization redshifts. Although the 
L r -R e and the Faber- Jackson relation do not pro- 
vide independent tests of the model (see Sect. 3.3) 
it is reassuring that two different data sets con- 
sistently confirm that it yields the correct virial 
mass-to-light ratio. 

The evolution of the stellar component is tied 
to that of the central black hole, and therefore 
the Mbh _ c relation is also shaped by the effect of 
feedback, primarily from the active nucleus itself 



in the more massive systems, and from supernovae 
in smaller objects. Again, observations are well 
reproduced (see Fig. 4), the dispersion around the 
best-fit relation is expected to be largely due to the 
different virialization redshifts, and a steepening 
of the Mbh~ct is expected for low BH masses, if 
these are mostly due to accretion interrelated with 
the bulge formation. However, small BH masses 
may be substantially increased by later accretion, 
occurring when the host galaxy is disk-dominated. 
Moreover, according to some analyses, seed BH 
masses may be ~ 10 5 M Q , and the low-tr portion 
of the diagram is testing more the distribution of 
seed masses than the accretion history. 
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We have further pointed out that the Mbh - 
M sp h relationship is essentially independent of z v - lr 
(see Fig. 5); its intrinsic scatter should therefore 
be minimum, although this docs not necessarily 
translate in a low observed scatter due to the large 
uncertainties on the M sp h estimates (compared to 
those on a). 

A synoptic view of the effect of feedback on the 
luminosity-er relation, on the local r*-band lumi- 
nosity function, and on the Mbh - c relation, based 
on the ABC model, is provided by Fig. 6. If we 
switch off the supernova feedback (central panel 
row), we get a larger luminosity at fixed virial 
velocity Kir- The increase is of ~ 1 mag. for 
low luminosity/mass galaxies (M r * ~ —20), where 
the stellar feedback is dominating over the QSO 
feedback, and of 0.5 mag. at high luminosities 
(M r . ~ -23). Turning off the QSO feedback (bot- 
tom panel row) goes in the same direction, but 
now the effect is larger for larger galaxies. To fit 
the observed Faber- Jackson relation with no stel- 
lar or QSO feedback, we have to shift by 0.1 dex 
the cr/Kir ratio, but this is inconsistent with the 
observed VDF. 

As shown by the central panel column of Fig. 6, 
the shift to higher luminosities occurring when the 
feedback is switched off affects only weakly the low 
luminosity portion of the luminosity function, be- 
cause of its flat slope. On the other hand, the high 
luminosity tail is very sensitive to it, and particu- 
larly to the feedback from the active nuclei. 

On the whole, the local luminosity function of 
galaxy spheroids and the observed correlations be- 
tween their properties provide clear evidence that 
the feedback both from supernovae and from ac- 
tive nuclei plays a key role in the evolution of these 
sources, and yield rather stringent constraints on 
the parameters controlling the coupling of the en- 
ergy injected into the interstellar medium. 
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